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Abstract 

Inclusive-jet cross sections have been measured in the reaction ep — >■ e + jet + X 
for photon virtuality < 1 GeV^ and 7p centre-of-mass energies in the region 
142 < W^p < 293 GeV with the ZEUS detector at HERA using an integrated 
luminosity of 300 pb~^. Jets were identified using the kr, anti-/cT or SIScone 
jet algorithms in the laboratory frame. Single-differential cross sections are pre- 
sented as functions of the jet transverse energy, -E;^*, and pseudorapidity, r/j'^*, 
for jets with i?;^* > 17 GeV and — 1 < 77^°* < 2.5. In addition, measurements of 
double-differential inclusive-jet cross sections are presented as functions of E^r^^ 
in different regions of rf'^^. Next-to- leading-order QCD calculations give a good 
description of the measurements, except for jets with low E}^^ and high 77-''^*. The 
influence of non-perturbative effects not related to hadronisation was studied. 
Measurements of the ratios of cross sections using different jet algorithms are 
also presented; the measured ratios are well described by calculations including 
up to 0{al) terms. Values of as(M^) were extracted from the measurements 
and the energy-scale dependence of the coupling was determined. The value 
of as{Mz) extracted from the measurements based on the kx jet algorithm is 
as{Mz) = 0.1206lo;Q022(6xp.)!'!o oo35(th-)! results from the anti-fc^ and SIS- 
cone algorithms are compatible with this value and have a similar precision. 
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1 Introduction 



The study of jet production in ep collisions at HERA has been well established as a testing 
ground of perturbative QCD (pQCD). Jet cross sections provided precise determinations 
of the strong coupling constant, as, and its scale dependence. The jet observables used 
to test pQCD included inclusive-jet [M], dijet [DIllEHS] and multijet [6l[71[l0Hl2] cross 
sections in neutral current (NC) deep inelastic ep scattering (DIS), inclusive-jet [T^IH], 
dijet [T5ti20] and multijet [2T1|22] cross sections in photoproduction and the internal struc- 
ture of jets in NC and charged current PM27] DIS. These studies also demonstrated 
that the cluster algorithm [28] in the longitudinally invariant inclusive mode [29] re- 
sults in the smallest uncertainties in the reconstruction of jets in ep collisions. Jet cross 
sections in NC DIS [2] and photoproduction [TH] were used by ZEUS [SU] as input in a 
QCD analysis to extract the parton distribution functions (PDFs) of the proton; these 
data helped to constrain the gluon density at medium- to high- a; values, where x is the 
fraction of the proton momentum carried by the gluon. 

The kx algorithm is well suited for ep collisions and yields infrared- and collinear-safe cross 
sections at any order of pQCD. However, it might not be best suited to reconstruct jets in 
hadron-hadron collisions, such as those at the LHC. In order to optimise the reconstruction 
of jet observables in such environments, new infrared- and collinear-safe jet algorithms 
were recently developed, namely the anti-/cT [21], a recombination-type jet algorithm, and 
the "Seedless Infrared-Safe" cone (SIScone) [32] algorithms. Measurements of jet cross 
sections in NC DIS using these algorithms were recently published [33] and constituted 
the first measurements with these new jet algorithms. The results tested the performance 
of these jet algorithms with data in a well understood hadron-induced reaction and it was 
shown that pQCD calculations with up to four partons in the final state provide a good 
description of the differences between jet algorithms. 

Measurements of inclusive-jet cross sections in photoproduction are presented in this 
paper. Two types of QCD processes contribute to jet production in photoproduction; 
at leading order they can be separated into [3ll|35] the direct process, in which the 
photon interacts directly with a parton in the proton, and the resolved process, in which 
the photon acts as a source of partons, one of which interacts with a parton in the 
proton. Due to the presence of the resolved processes, the analysis of jet cross sections 
in photoproduction with different jet algorithms provides a test of their performance in a 
reaction closer to hadron-hadron interactions than NC DIS. 

In this paper, single-differential inclusive-jet cross sections as functions of the jet trans- 
verse energy, -E^*, and pseudorapidity, rj^^^, are presented based on the /cy, anti-fc^ and 
SIScone jet algorithms. The results based on the anti-Zcr and SIScone jet algorithms are 
compared to the measurements based on the kr via the ratios of cross sections. In ad- 
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dition, measurements of cross sections are also presented as functions of E';^ in different 
regions of tf^^, which have the potential to constrain further the gluon density at high 
X. Next-to-leading-order (NLO) QCD calculations using recent parameterisations of the 
proton and photon PDFs are compared to the measurements. A determination of as{Mz) 
as well as of its energy-scale dependence are also presented. The analyses presented here 
are based on a data sample with a more than three-fold increase in statistics with respect 
to the previous study [13]. 



2 Experimental set-up 

A detailed description of the ZEUS detector can be found elsewhere [361 EZ]- A brief 
outline of the components most relevant for this analysis is given below. 

Charged particles were tracked in the central tracking detector (CTD) [SB], which operated 
in a magnetic field of 1.43 T provided by a thin superconducting solenoid. The CTD 
consisted of 72 cylindrical drift-chamber layers, organised in nine superlayers covering the 
polar-angl^ region 15° < 9 < 164°. The CTD was complemented by a silicon microvertex 
detector (MVD) [39], consisting of three active layers in the barrel and four disks in the 
forward region. For CTD-MVD tracks that pass through all nine CTD superlayers, the 
momentum resolution was o"(pt)/pt = 0.0029pT © 0.0081 © 0.0012/pt, with px in GeV. 

The high-resolution uranium-scintillator calorimeter (CAL) [40] consisted of three parts: 
the forward (FCAL), the barrel (BCAL) and the rear (RCAL) calorimeters. Each part was 
subdivided transversely into towers and longitudinally into one electromagnetic section 
(EMC) and either one (in RCAL) or two (in BCAL and FCAL) hadronic sections (HAC). 
The smallest subdivision of the calorimeter was called a cell. The CAL energy resolutions, 
as measured under test-beam conditions, were a{E)/E = 0.18/ \^ for electrons and 
a{E)/E = 0.35//E for hadrons, with E in GeV. 

The luminosity was measured using the Bethe-Heitler reaction ep — e'-fp by a luminos- 
ity detector which consisted of a lead-scintillator calorimeter [H] and an independent 
magnetic spectrometer [12]. The fractional uncertainty on the measured luminosity was 
1.8%. 



The ZEUS coordinate system was a right-handed Cartesian system, with the Z axis pointing in the 
proton beam direction, referred to as the "forward direction" , and the X axis pointing towards the 
centre of HERA. The coordinate origin was at the nominal interaction point. 
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3 Data selection 



The data were collected during the running period 2005-2007, when HERA operated with 
protons of energy Ep = 920 GeV and electrons or positron^ of energy Ef, = 27.5 GeV, at 
an ep centre-of-mass energy of y/s = 318 GeV, and correspond to an integrated luminosity 
of 299.9 ±5.4 ph~\ 

A three-level trigger system was used to select events online [371113]. At the first level, 
events were triggered by a coincidence of a regional or transverse energy sum in the GAL 
and at least one track from the interaction point measured in the GTD. At the second 
level, a total transverse energy of at least 12 GeV, excluding the energy in the eight 
GAL towers immediately surrounding the forward beampipe, was required, and cuts on 
GAL energies and timing were used to suppress events caused by interactions between the 
proton beam and residual gas in the beampipe. At the third level, two different methods 
were applied to select the events. The first method selected events with a total transverse 
energy of at least 25 GeV, excluding the energy in the eight GAL towers immediately 
surrounding the forward beampipe. For the second method, a jet algorithm was applied 
to the GAL cells and jets were reconstructed using the energies and positions of these 
cells; events with at least one jet of Et > 10 GeV and r] < 2.5 were accepted. Additional 
requirements based on GAL energies, tracking and timing were used to suppress further 
the non-ep background. 

Events from collisions between quasi-real photons and protons were selected offline using 
similar criteria to those reported in the previous ZEUS publication [13]. The selection 
criteria applied were: 

• a reconstructed event vertex along the Z axis within 35 cm of the nominal interaction 
point was required; 

• cuts based on tracking information were applied to remove the contamination from 
beam-gas interactions, cosmic-ray showers and beam-halo muons; 

• charged current DIS events were rejected by requiring the total missing transverse 
momentum, p^^^^, to be small compared to the total transverse energy, E^^, i.e. 
pf^'^/^E^ < 2 VGeV; 

• any NG DIS event with an identified scattered-electron candidate in the GAL was 

rejected; 

• the events were restricted to 7p centre-of-mass energies in the region 142 < W^p < 293 GeV, 
where W^p = ^/sy; y is the inelasticity and was estimated as yjB = {E — pz)/2Ee, 

^ In the following, the term electron will refer to both the electron and positron, unless otherwise stated. 
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where E is the total energy measured in the CAL and pz is the longitudinal component 
of the total momentum. 

After these selection criteria were applied, the contamination from beam-gas interactions, 
cosmic-ray showers and beam-halo muons was found to be negligible. The remaining 
background from NC DIS events was estimated by Monte Carlo (MC) techniques to be 
around 1% and was neglected. The contamination from charged current DIS events was 
found to be even smaller. The selected sample consisted of events from ep interactions 
with < 1 GeV^, where is the virtuality of the exchanged photon, and a median 
Q'^ ^ 10~^ GeV^, estimated using MC techniques. 

4 Jet search 

In photoproduction, jets are usually defined using the transverse-energy fiow in the 
pseudorapidity-azimuth [rj — (p) plane of the laboratory frame. The procedure to recon- 
struct jets with the algorithm from an initial list of objects (e.g. final-state partons, 
final-state hadrons or energy deposits in the calorimeter) is described below in some de- 
tail. In the following discussion, denotes the transverse energy, the pseudorapidity 
and 0* the azimuthal angle of object i. For each pair of objects, the quantity 

d,, = minmr, ■ [{rj^ - r^^f + (0^ - (1) 

is calculated, where R is the jet radius. For each individual object, the distance to the 
beam, di = (E^)"^, is also calculated. If, of all the values {dij, di}, dki is the smallest, then 
objects k and / are combined into a single new object. If, however, dk is the smallest, 
then object k is considered a jet and removed from the sample. The procedure is repeated 
until all objects are assigned to jets. 

The anti-/cr algorithm is identical to the kx except for a modified distance measure, 

d,, = min((E^)-2, . [(rj^^ - r^^f + (0^^ - ^p^fyR^ (2) 

and the distance to the beam, which is defined as di = 

The SIScone algorithm consists of two steps. First, for a given set of initial objects, all 
stable cones are identified; cones are classified as stable by the coincidence of the cone 
axis with that defined by the total momentum of the objects contained in the given cone 
of radius R in the rj — cj) plane. In this procedure, no seed is used. Stable cones are 
then discarded if their transverse momentum is below a given threshold, pt,mm- For each 
selected stable cone, the scalar sum of the transverse momentum of the objects associated 
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to it, pt, is defined. Second, overlapping cones are identified and subsequently split or 
merged according to the following procedure. Two cones are merged if the scalar sum 
of the transverse momentum of the objects shared by the two cones exceeds a certain 
fraction / of the lowest-pt cone; otherwise, two different cones are considered and the 
common objects are assigned to the nearest cone. 

For the measurements presented in this paper, the jet radius R was set to unity and 
the jet variables were defined according to the Snowmass convention [H] for all three jet 
algorithms. In the application of the SIScone algorithm, the fraction / was set to 0.75 

and pt, min was set to zero. 

The kx, anti-/cT and SIScone jet algorithm^ were used to reconstruct jets in the hadronic 
final state from the energy deposits in the CAL cells. The jets reconstructed from the 
CAL cell energies are called calorimetric jets and the variables associated with them are 
denoted by -Ey ^ai' ''7cai '/'cai- Three samples of events were selected for further analysis, 
one for each jet algorithm, which contain at least one jet satisfying -E'^cai ^ GeV and 
-1 < < 2.5. 

5 Monte Carlo simulations 

Samples of events were generated to determine the response of the detector to jets of 
hadrons and the correction factors necessary to obtain the hadron-level jet cross sections. 
In addition, these samples were used to estimate hadronisation corrections to the NLO 
calculations (see Section [H]). 

The MC programs Pythia 6.146 |16] and Herwig 6.504 |17] were used to generate 
resolved and direct photoproduction events. In both generators, the partonic processes 
are simulated using leading-order matrix elements, with the inclusion of initial- and final- 
state parton showers. Fragmentation into hadrons was performed using the Lund string 
model IH] , as implemented in Jetset [19] in the case of Pythia, and a cluster model [50] 
in the case of Herwig. The CTEQ4M [51] (GRV-HO [S2]) sets were used for the proton 
(photon) PDFs. Samples of Pythia including multiparton interactions [53] (Pythia- 
MI) with a minimum transverse momentum for the secondary scatter, p^^jn, of 1, 1.5 or 
2 GeV were used to simulate contributions from non-perturbative effects not related to 
hadronisation (NP), such as the underlying event. All the samples of generated events 
were passed through the Geant 3.21-based [SI] ZEUS detector- and trigger-simulation 
programs [37]. They were reconstructed and analysed by the same program chain as the 
data. 

^ The Fastjet 2.4.1 [JS] package was used. 
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The jet search was performed on the MC events using the energy measured in the CAL 
cells as described in Section HJ In addition, the three jet algorithms were also applied 
to the final-state particles (hadron level) and partons (parton level). The hadron level is 
defined by those hadrons with lifetime r > 10 ps and the parton level is defined as those 
partons present after the parton-shower procedure. 

6 Transverse-energy and acceptance corrections 

The comparison of the reconstructed jet variables for the hadronic and the calorimetric jets 
in MC-simulated events showed that no correction was needed for the jet pseudorapidity 
and azimuth. However, -E^*^! underestimates the corresponding hadronic-jet transverse 
energy by ~ 14% with an r.m.s. of ~ 10%. This underestimation is mainly due to the 
energy lost by the particles in the inactive material in front of the CAL. The transverse- 
energy corrections to calorimetric jets, as functions of 77^°*^ and E^rj^^^^ and averaged over 
(jy'^^i, were determined using the MC events. Further corrections to the jet transverse 
energy were applied to the data to account for differences in the jet energy scale between 
data and MC simulations; the method presented previously pTtlSB] , which relies on a good 
understanding of the performance of the track reconstruction, was used to calibrate the 
absolute energy scale of the jets down to ±1%. This calibration was cross checked by 
means of the transverse-momentum balance in the CAL between the electron candidate 
and the jet in single-jet NC DIS events. 

Henceforth, jet variables without subscript refer to the corrected values. After all these 
corrections to the jet transverse energy, events with at least one jet satisfying E'^* > 
17 GeV and —1 < 77-''^* < 2.5 were retained. The number of events and jets in the final 
data samples are shown in Tabled] for each jet algorithm. No events with more than four 
jets were found in these samples. 

The E^r^^ and 7]^°^ distributions in the data were corrected for detector effects using bin- 
by-bin acceptance correction factors determined using the MC samples. These correction 
factors take into account the efficiency of the trigger, the selection criteria and the purity 
and efficiency of the jet reconstruction. For this approach to be valid, the uncorrected dis- 
tributions of the data must be adequately described by the MC simulations at the detector 
level. This condition was satisfied by both the Pythia and Herwig MC samples. The 
average between the acceptance correction factors obtained from Pythia and Herwig 
was used to correct the data to the hadron level. The deviations in the results obtained 
by using either Pythia or Herwig to correct the data from their average were taken 
to represent systematic uncertainties of the effect of the QCD-cascade and hadronisation 
models in the corrections (see Section [7]). The acceptance correction factors differed from 
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unity by typically less than 20%. 



7 Experimental uncertainties 

The following sources of systematic uncertainty were considered for the measured cross 
sections: 

• the differences in the results obtained by using either Pythia or Herwig to correct 
the data for detector effects. The resulting uncertainty was typically below ±4%; 

• the effect of the CAL energy-scale uncertainty on W^^p was estimated by varying yjB 
by ±1% in simulated events. The uncertainty in the cross sections was below ±1% at 
low -E;^*, increasing to ~ ±3% at high -E^*; 

• the effect of the uncertainty on the parameterisations of the proton and photon PDFs 
was estimated by using alternative sets of PDFs (MRST(c-g) for the proton and 
AFG-HO [57j for the photon) in the MC simulation to compute the acceptance cor- 
rection factors. The variation of the cross sections was typically smaller than ±1%; 

• the uncertainty in the cross sections due to that in the simulation of the trigger was 
found to be negligible. 

All the above systematic uncertainties were added in quadrature and are shown in the 
figures as error bars. The resulting systematic uncertainty in the cross sections based on 
the three jet algorithms was similar and typically below ±5%. 

The absolute energy scale of the calorimetric jets in simulated events was varied by its 
uncertainty of ±1% (see Section E]); the effect of this variation on the inclusive-jet cross 
sections was typically =f5% at low -E^*, increasing up to ^10% at high E}^\ This uncer- 
tainty is fully correlated between measurements in different bins and is shown separately 
as a shaded band in the figures. In addition, there was an overall normalisation uncer- 
tainty of ±1.8% from the luminosity determination, which is not included in the figures 
and tables. 

8 Next-to-leading-order QCD calculations 

The NLO QCD ((9(a^)) calculations used in the analysis presented here were computed 
using the program by Klasen, Kleinwort and Kramer [58]. The calculations use the phase- 
space-slicing method [52] with an invariant-mass cut to isolate the singular regions of the 
phase space. The number of flavours was set to five and the renormalisation {fin) and 
factorisation (fxp) scales were chosen to he = I^f = 1^ = E^t^. The strong coupling 
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constant was calculated at two loops with = 226 MeV, corresponding to as{Mz) = 
0.118. The calculations were performed using the ZEUS-S [HO] parameterisations of the 
proton PDFs and the GRV-HO sets for the photon PDFs as defaullQ. The three jet 
algorithms were applied to the partons in the events generated by this program to compute 
the jet cross-section predictions. At 0{a'j.), the parton-level predictions from the 
and anti-/cT jet algorithms are identical. The predictions from the kr and SIScone jet 
algorithms start to differ at this order. 

Since the measurements refer to jets of hadrons, whereas the NLO QCD calculations refer 
to jets of partons, the predictions were corrected to the hadron level using the MC models. 
The multiplicative correction factor, Chad; was defined as the ratio of the cross section 
for jets of hadrons over that for jets of partons, estimated by using the MC programs 
described in Section O The mean of the ratios obtained with Pythia and Herwig was 
taken as the value of Chad- Details on the values of Chad are presented in Section [91 

The following sources of uncertainty in the theoretical predictions were considered: 

• the uncertainty on the NLO QCD calculations due to that on the value of as{Mz) 
used was estimated by repeating the calculations using two additional ZEUS-S sets 
of proton PDFs, for which different values of as{Mz) were assumed in the fits. The 
difference between these calculations was used to determine the uncertainty due to 
that on the value of as{Mz) [62] : 

• the uncertainty on the NLO QCD calculations due to terms beyond NLO was estimated 
by repeating the calculations using values of /i scaled by factors 0.5 and 2; 

• the uncertainty from the modelling of the QCD cascade and hadronisation effects was 
assumed to be half the difference between the hadronisation corrections obtained using 
the Pythia and Herwig models; 

• the uncertainty on the NLO QCD calculations due to those on the proton PDFs was 
estimated by repeating the calculations using 22 additional sets from the ZEUS-S error 
analysis with as{Mz) fixed to the central value; this error analysis takes into account 
the statistical and correlated experimental uncertainties of each data set used in the 
determination of the proton PDFs; 

• the uncertainty on the NLO QCD calculations due to those on the photon PDFs was 
estimated by using alternative sets of parameterisations, AFG04 [63] and CJK [6^ . 

The total theoretical uncertainty was obtained by adding in quadrature the individual 
uncertainties listed above. Figure [1] shows an overview of the relative theoretical uncer- 
tainties for the inclusive-jet cross sections in the kinematic region of the measurements as 

The Lhapdf 5.7.1 [HI] package was used. 
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functions of rj^^^ and E^r^^ and for each jet algorithm separately. The uncertainty due to 
higher orders is somewhat larger for the SIScone than for the /ct and anti-Zc^ algorithms, 
whereas the other uncertainties are very similar for the three jet algorithms. The uncer- 
tainty coming from the terms beyond NLO is dominant in all cases. At high the 
proton PDF uncertainty is of the same order (slightly smaller) as that arising from terms 
beyond NLO for the fcy and anti-fcy (SIScone) algorithms. The uncertainty arising from 
the photon PDFs at high 77-^''* becomes comparable to that coming from higher orders. 
The uncertainties from the value of as{Mz) and hadronisation corrections are small. 

The samples of Pythia-MI described in Section [5] were used to estimate the contribution 
from non-perturbative effects not related to hadronisation. Such effects were computed 
as ratios of the cross section for jets of hadrons in the Pythia-MI samples over that 
for the samples of Pythia; these ratios are called Cnp. The values of Cnp depend 
strongly on Px'min^ ^^e minimum transverse momentum for the secondary scatter, since 
the smaller the Pr,min value is set, the larger the phase space available for production 
of secondary interactions and hence the higher the jet rate. Another feature of such 
secondary interactions is that, due to the particular kinematics of HERA, the products 
of these additional interactions are expected to be boosted towards the proton direction. 
Several predictions including these non-perturbative effects, denoted as NLO(g>NP, were 
computed by applying the factors Cnp, using the Prmin — 1' ^-^ 2 GeV Pythia-MI 
samples, to the NLO QCD calculations after hadronisation corrections. 

9 Results 

Single- and double-differential inclusive-jet cross sections were measured in the kinematic 
region given by < 1 GeV^ and 142 < W:yp < 293 GeV. These cross sections include 
every jet of hadrons with i?;^* > 17 GeV and —1 < rj^^^ < 2.5 in each event. The jets 
were reconstructed using either the kx, the anti-/cT or the SIScone jet algorithms. The x 
region covered by the measurements was determined to be 3 ■ 10~^ < x < 0.95. 

9.1 Single-differential cross sections 

The measurements of the single-differential cross sections based on the /ct jet algorithm 
as functions of -E;^* and rj^^^ are presented in Figs. [2] and [3] and Table [21 In these and 
the subsequent figures, each data point is plotted at the weighted mean of each bin. The 
measured da/dEi^^ falls by over four orders of magnitude in the measured range. The 
measured da/drf'^^ displays a maximum around rf^^ ~ 1. 
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The NLO QCD predictions are compared to the measurements in these figures. The 
calculation reproduces the measured da/dE}^^ well. The measured da/drf^^^ is well de- 
scribed for rf^^ < 2; however, an excess of data with respect to the theory is observed for 
larger t]-"^* values. Such discrepancies have already been observed in previous studies of 
jet photoproduction [^DIES]- In the study of dijet production [2U], the discrepancies were 
interpreted as an inadequacy of the parameterisations of the photon PDFs, which had 
been extracted from e^e~ data at lower scales. The studies of multijet production [22] 
showed the need to include non-perturbative effects not related to hadronisation in the 
pQCD calculations to describe the data. 

The infiuence of non-perturbative effects not related to hadronisation in the predictions 
was investigated by using the NLO®NP QCD calculations (see Section [S]). The compari- 
son of these calculations to the data is shown in Fig. HJ It is observed that the NLO®NP 
QCD calculations predict a larger jet rate at low E^^^ and high Tf'^^, in the region where 
the NLO QCD prediction fails to describe the data. The NLO®NP QCD prediction with 
PT,min = 1-5 GeV is closest to the data. These observations indicate the possible presence 
of effects such as the underlying event in the data, which are not included in the NLO 
QCD calculation. These non-perturbative contributions are expected to be unrelated 
to the hard scattering and approximately constant with the scale of the interaction, so 
that the ratio of this non-perturbative contribution to the jet transverse energy becomes 
smaller as E^^^ increases, as seen in Fig. H^. This is supported by the good description of 
the data by the NLO QCD calculation for E-!^'^ > 21 GeV (Fig.S^) and by the inclusive-jet 
cross section function of r/J'^* for Ei^' > 21 GeV (see Fig. [5] and Table [3]): the NLO 
QCD calculation gives a good description of the data in the whole measured range; in 
particular, discrepancies between data and theory are no longer observed at high values 
of Tyj*^*. In addition, the differences between the NL0(8)NP predictions with different Prmm 
values become smaller, as seen in Fig. 03. 

The influence of the poorly constrained photon PDFs on the predictions was investigated 
by comparing calculations based on different PDF sets to the data. Figure [6] shows the 
measurements together with the NLO QCD predictions using alternatively the AFG04 
and CJK sets of photon PDFs, together with the predictions based on the GRV-HO set. 
Some differences are observed between the three predictions, especially at low i?;^* and 
high 77^'^*. In particular, the predictions based on AFG04 (CJK) are lower (higher) than 
those based on GRV-HO. 

The influence of the proton PDFs on the predictions was investigated by comparing cal- 
culations based on different PDF sets to the data. Figure [7] shows the measurements 
together with the NLO QCD predictions using alternatively the MSTW08 [HS] and HER- 
APDF1.5 [66] sets of proton PDFs, together with the predictions based on the ZEUS-S 
set. The prediction based on HERAPDF1.5 is lower than that based on the ZEUS-S set 
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in most of the phase space, whereas the MSTW08 prediction is higher at high E^r^ . This 
region of phase space is not well constrained since the main contribution comes from the 
high-x gluon density in the proton. 

In summary, the measurements of inclusive-jet cross sections in photoproduction have 
the potential to constrain the proton and the photon PDFs. To study in more detail 
the sensitivity of the inclusive-jet cross sections to the proton and photon PDFs and find 
the regions of phase space in which the data can add information to constrain further 
these PDFs, double-differential cross sections were measured and are presented in the 
next section. 

9.2 Double-differential cross sections 

The measurements of the inclusive-jet cross sections based on the kj- jet algorithm as 
functions of -E^* in different regions of rj^^^ are presented in Fig. |8] and Tables H] and O 
The measured cross sections exhibit a steep fall-off within the E';^* range considered. The 
E^r^^ dependence of the cross section becomes less steep as r]^^^ increases. 

The NLO QCD predictions are compared to the measurements in Fig. [HI They give a 
good description of the data, except at low i?;^* and high t]^^^. Figure IH] shows the relative 
difference of the measured differential cross sections to the NLO QCD calculations. The 
data are well described by the predictions for —1 < rj^^^ < 2 in the whole E^r^^ range 
measured. For the region 2 < 77^'^* < 2.5, where it is observed that non-perturbative 
effects not related to hadronisation (see Fig. [TOj) might contribute significantly, the data 
are well described only for E";^* > 21 GeV. 

Figures [TT] and [T^] show the comparison between the measured cross sections and the 
predictions based on different photon and proton PDFs, respectively. As discussed above, 
differences at low E";^* and high t]^^^ are observed between the predictions based on GRV- 
HO, AFG04 and CJK. The latter gives predictions closest to the data, especially in 
the region 2 < rf^^ < 2.5. The largest differences between the predictions based on 
MSTW08 and ZEUS-S are observed at high E^r^^ for r]^^^ > 1. The predictions based on 
HERAPDF1.5 are lower than those based on ZEUS-S in most of the phase-space region. 

As discussed in Section [8], the theoretical uncertainties are dominated by the contribution 
from higher orders. This uncertainty decreases as increases. The contribution from 
the proton PDF uncertainty is significant and approximately constant for i?^* > 30 GeV; 
at high E^r^^ values, the proton PDF uncertainty is of the same order as that coming from 
higher orders. In these regions, in which the gluon-induced contribution is still substantial 
and the possible presence of non-perturbative effects is expected to be minimised, the 
data have the potential to constrain the gluon density in the proton. The uncertainty 
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coming from the photon PDFs is largest at low E^r^^ and high rj^^^ and approximately of 
the same order as that coming from higher-order terms. Therefore, these high-precision 
measurements also have the potential to constrain the photon PDFs in these regions of 
phase space. 

9.3 Single-differential cross sections based on different jet algo- 
rithms 

The measurements of the inclusive-jet cross sections based on the anti-/cT and SIScone jet 
algorithms as functions of E^r^^ and 77-^°* are presented in Fig. [131 together with those based 
on the kx algorithm. The measured cross sections are also given in Tables M and [71 The 
measured da /dE^^^ cross sections exhibit a steep fall-off of over four orders of magnitude 
in the E'y* measured range. The measured da/drf^^ cross sections display a maximum 
around 77-^''* ~ 1. The measured cross sections using the three jet algorithms have a similar 
shape, normalisation and precision. 

The NLO QCD predictions are compared to the data in Fig. [121 The hadronisation 
correction factors applied to the calculations and their uncertainties are also shown. It 
is seen that the hadronisation correction factors are closest to (farthest from) unity for 
the kx (SIScone) jet algorithm (see also Tables [SI and [7|) . The ratios of the measured 
cross sections to the NLO QCD calculations are shown in Fig. [TH separately for each jet 
algorithm. The measured cross sections are well reproduced by the calculations, except 
at high 77-'°*. 

The ratios of the cross-sections anti-Zcr/^T, SIScone//cr and anti-/cr/SIScone were studied 
to compare the jet algorithms in more detail. These ratios allow, in particular, a stringent 
test of the description of the differences between jet algorithms in terms of parton radi- 
ation due to the partial cancellation of experimental and theoretical uncertainties. The 
measured ratios are shown in Fig. [TJl In these ratios, the statistical correlations among 
the event samples as well as those among the jets in the same event were taken into 
account in the estimation of the statistical uncertainties; most of the systematic uncer- 
tainties, including that due to the jet energy scale, cancel out. The measurements show 
that the cross sections based on the anti-Zcy algorithm are similar in shape to those based 
on the kT algorithm but ~ 6% lower and that the cross sections based on the SIScone 
have a slightly different shape than those based on the kx and anti-/cr algorithms. 

The pQCD predictions including up to 0{a1) terms for the ratios of the cross sections 
are also shown in Fig. [13 In the estimation of the total theoretical uncertainty of the 
predicted ratios, all the theoretical contributions were assumed to be correlated except 
those due to terms beyond 0{a1) and to the modelling of the QCD cascade and hadro- 
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nisation. Figure [15] also includes the ratio of the hadronisation correction factors applied 
to the 0(q!^) predictions. The predictions for the anti-Zc^ and kx algorithms are iden- 
tical at 0{a1) and, thus, the predicted anti-Zcr/^T ratio coincides with the ratio of the 
hadronisation corrections for both algorithms. The predictions for the kr and anti-/cr 
algorithms are expected to start to differ at 0{a^) and, conservatively, the uncertainties 
due to higher-order terms were assumed to be uncorrelated in the ratio; otherwise, a co- 
herent variation of the renormalisation scale in the ratio would yield an unrealistic zero 
contribution. In the case of the SlScone/Zcr and anti-Zcr/SIScone, the C(q;^) predictions 
are different and the ratios differ significantly from the ratios of the hadronisation correc- 
tion factors. The measured ratios are well described by the calculations including terms 
up to 0{a^) within the experimental and theoretical uncertainties, which are dominated 
by the uncertainty due to higher orders. 

The differences in the influence of non-perturbative effects not related to hadronisation 
for the different jet algorithms were investigated by comparing the ratios of the NLO®NP 
predictions with Prmin — ^-^ GeV, which are also shown in Fig. [151 For the anti-/cT/^r 
ratio, the NLO and the NLO®NP predictions are very similar, which indicates that the 
non-perturbative effects affect the jets in the same way. However, in the SlScone/Zc^ and 
anti-Zcr/SIScone ratios, the ratios of the NLO(g)NP predictions differ from the ratios of 
the NLO calculations at low E^r^^ and high t]^^^; these differences, which are located in the 
regions of phase space where the NLO calculations fail to describe the data, show that 
there is a dependence on the details of the jet reconstruction concerning non-perturbative 
effects not related to hadronisation. These differences are at most of the same size as the 
theoretical uncertainties. 

In summary, it is concluded that the measured inclusive-jet cross sections based on the 
/ct, anti-fcr and SIScone jet algorithms are well described by the NLO QCD calculations, 
except at high rf'^^. The data and the calculations for the three jet algorithms have a 
similar experimental and theoretical precision. Furthermore, the measured ratios are well 
described by the predictions including up to O(a^) terms, demonstrating the ability of 
the pQCD calculations with up to three partons in the final state to account adequately 
for the differences between the SIScone and the kx or anti-/cT jet algorithms. 

9.4 Determination of a.s{Mz) 

The measured single-differential cross-sections da/dE^r^^ based on the three jet algorithms 
were used to determine values of as{Mz) using the method presented previously [2]. The 
NLO QCD calculations were performed using five different sets of the ZEUS-S proton 
PDFs which were determined from global fits assuming different values of as{Mz), namely 
as{Mz) = 0.115, 0.117, 0.119, 0.121 and 0.123. The GRV-HO set was used as default for 
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the photon PDFs. The value of as{Mz) used in each calculation was that associated with 
the corresponding set of proton PDFs. 

The as{Mz) dependence of the predicted cross sections in each bin i of -E;^* was parame- 
terised according to 

da/dE^^\as{Mz))] = C{a,{Mz) + ^^^(m^), 

J i 

where C\ and C'2 were determined from a fit to the NLO QCD calculations. The value 
of as{Mz) was determined by a fit to the measured da/dE^^^ values. In the fitting 
procedure, the running of as predicted by QCD was assumed. Only the measurements 
for E^^ > 21 GeV were used in the fit to minimise the effects of a possible non-perturbative 
contribution in addition to that of hadronisation and the uncertainty coming from higher 
orders. In addition, the fit was restricted to i?^* < 71 GeV because of the relatively large 
uncertainty coming from the proton PDFs for higher i?;^* values. 

The experimental uncertainties on the extracted values of as{Mz) were evaluated by 
repeating the analysis for each systematic check presented in Section [71 The overall 
normalisation uncertainty from the luminosity determination was also included. The 
largest contribution comes from the uncertainty in the absolute jet energy scale. The 
theoretical uncertainties were evaluated as described in Section [SJ The largest contribution 
arises from the terms beyond NLO, which was estimated by using the method of Jones 
et al. [HZ]. The uncertainty due to the photon PDFs is of the same order as that arising 
from higher orders. The uncertainty due to the proton PDFs and that arising from 
the hadronisation effects were also estimated. All uncertainties are listed separately in 
Table El 

As a cross-check, as{Mz) was determined by using NLO QCD calculations based on the 
CTEQ6.1 [68] (MSTW08) sets of proton PDFs. The values obtained are consistent within 
1.0 (1.0), 0.9 (0.9) and 0.6 (0.8)% with those based on ZEUS-S for the kr, anti-^T and SIS- 
cone determinations, respectively. The uncertainty arising from the proton PDFs was esti- 
mated to be ±1.3 (0.9, 0.5)% for the kr and anti-fcr and ±1.2 (0.8, 0.4)% for the SIScone 
determinations using the results of the CTEQ6.1 (MSTW08nlo90cl, MSTW08nlo68cl) 
analysis. 

The values of as{Mz) obtained from the measured da /dE^^^ are 



as{Mz)K = 0.1206 1°:°°^^ (exp.) 1™ (th.), 
as{Mz)U-,, = 0.1198 t'ofoll (exp.) Tootl (th-), 
«.(M^)|sisco„e = 0.1196 l°:°°^2(exp.)+™(th.). 
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The value of as{Mz) determined from the anti-Z^T (SIScone) measurements is consistent 
with that obtained from the kj- analysis within —0.7 (—0.8)%. The total uncertainty from 
each determination, as shown in Table [H indicates that the performance of the three jet 
algorithms is similar. These determinations are consistent with previous determinations in 
NC DIS [SH?] at HERA, with the results obtained in pp collisions pHI and have a precision 
comparable to those obtained in individual determinations from e~^e~ experiments |70j . 
These values are also consistent with the world average [H2], as well as with the HERA 
2004 average [71] and the HERA 2007 combined value [72]. Figure [16] shows the value of 
as{Mz) determined from the /c^-based analysis together with determinations from other 
experiments and the HERA and world averages. 

9.5 Energy-scale dependence of as 

The energy-scale dependence of was determined from a NLO QCD fit to the measured 
cross section. Values of were extracted at each mean value of measured E^^^ 
without assuming the running of Ug- The method employed was the same as described 
above, but parameterising the dependence of da/dE^^^ in terms of as{{E-!^^)) instead 
of as{Mz)., where (-E^*) is the average E^^^ of the data in each bin. 

The extracted values of function of E-!^^ are shown in Fig. [17] and Table [9] for the 

analysis based on the algorithm. The data demonstrate the running of over a large 
range in E';^* from a single experiment. The predicted running of the strong coupling [73] 
calculated at two loops is in good agreement with the data. The values of as functions 
of E^^^ determined from the anti-fcy and SIScone measurements are consistent with those 
obtained from the kx analysis and have a similar precision. 

10 Summary and conclusions 

Measurements of differential cross sections for inclusive-jet photoproduction at a centre-of- 
mass energy of 318 GeV using an integrated luminosity of 300 pb~^ collected by the ZEUS 
detector have been presented. The cross sections refer to jets of hadrons of E';^* > 17 GeV 
and —1 < 7]^^^ < 2.5 identified in the laboratory frame with the kx-, anti-fcy or SIScone jet 
algorithms with jet radius R= 1. The cross sections are given in the kinematic region of 
g2 < 1 GeV2 and 142 < W^p < 293 GeV. 

Measurements of single- differential cross sections were presented as functions of and 
j^jet rpj^g NLO QCD calculations provide a good description of the measured cross sections, 
except at high r/j^*. Non-perturbative effects not related to hadronisation and the influence 
of the photon PDFs were found to be most significant in this region. 
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Cross-section measurements were performed as functions of E^r^^ in different regions of r^j^*. 
The data are well described by the NLO QCD predictions, except for jets with low E^r^^ 
and high rf^^. These cross sections are sensitive to the parton densities in the proton and 
the photon in regions of phase space where the theoretical uncertainties are small. The 
precision measurements presented here are therefore of particular relevance for improving 
the determination of the PDFs in future QCD fits. 

A detailed comparison between the measurements for the three jet algorithms was per- 
formed. The measured cross sections for the three jet algorithms have a similar shape, 
normalisation and precision. The NLO QCD calculations of inclusive-jet cross sections 
and their uncertainties for the different jet algorithms were also compared: the predic- 
tions describe the data well, except at high 77^°*; the calculations based on the SIScone 
algorithm are somewhat less precise than those based on the /ct or anti-Zc-r due to a larger 
contribution from terms beyond NLO. The ratios of the cross sections based on the dif- 
ferent jet algorithms were also presented. The measured ratios are well reproduced by 
the C(tt^) predictions, demonstrating the ability of the pQCD calculations including up 
to three partons in the final state to account adequately for the details of the differences 
between the SIScone and the kr or anti-/cT jet algorithms. 

The measured cross sections were used to determine values of as{Mz). QCD fits to the 
cross-section da/dE^^^ for 21 < < 71 GeV based on the jet algorithm yielded 



a.(M^) = 0.1206 1™ (exp.) 1°°^^ (th.). 

This value is in good agreement with the world and HERA averages. The extracted values 
of CXg clS cl function of -Ey* are in good agreement with the predicted running of the strong 
coupling over a large range in E^^^ . 
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anti- /ct 


SIScone 


events 


483328 


444295 


438906 


jets 


613165 


572865 


566000 


one jet 


355691 


317509 


313519 


two jets 


125468 


125016 


123700 


three jets 


2138 


1756 


1667 


four jets 


31 


14 


20 



Table 1: Number of events and jets selected in data with E-!^^ > 17 GeV and 
— 1 < r/j''* < 2.5 in the kinematic region of Q'^ < I GeV^ and 142 < W^p < 293 
GeV for the kx, anti-kx and SIScone jet algorithms. The number of events with 
one, two, three and four jets are also listed. 
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E^r^^ bin 




da / dE^rp" 


'^stat 




<^ES 


Chad 


(GeV) 


(GeV) 


(pb/GeV) 










17 — 21 


IS 7 


295 84 


+0 52 


+6.37 
-6.36 


+ 11.35 
-11.93 


99 


21 - 25 


22.7 


95.86 


±0.28 


+1.56 
— 1.55 


+4.19 
-4.50 


0.99 


25 - 29 


26.7 


36.88 


±0.18 


+0.52 
-0.52 


+ 1.76 
-1.80 


0.98 


29-35 


31.4 


13.606 


±0.090 


+0.150 
-0.142 


+0.687 
-0.772 


0.98 


35-41 


37.5 


4.492 


±0.051 


+0.104 
-0.102 


+0.250 
-0.251 


0.98 


41-47 


43.6 


1.677 


±0.032 


+0.033 
-0.033 


+0.098 
-0.102 


0.98 


47-55 


50.2 


0.589 


±0.017 


+0.015 
-0.015 


+0.038 
-0.044 


0.98 


55-71 


60.3 


0.1216 


±0.0057 


+0.0033 
-0.0028 


+0.0086 
-0.0090 


0.97 


71-95 


77.2 


0.0121 


±0.0014 


+0.0010 
-0.0010 


+0.0012 
-0.0011 


1.00 





bin 






da/dv^"^ 
(pb) 


^stat 


'^syst 




Chad 


-0.75 - 


- -0 


50 


-0.59 


80.65 


±0.92 


+5.12 
-4.89 


+6.86 
-7.78 


0.82 


-0.50 - 


- -0 


25 


-0.36 


221.4 


±1.6 


+8.7 

-8.5 


+14.0 
-15.3 


0.93 


-0.25- 


-±o 


00 


-0.11 


404.4 


±2.2 


+10.0 
-9.9 


+21.7 
-23.1 


0.97 


±0.00 - 


- ±0 


25 


±0.13 


555.2 


±2.6 


+8.1 
-7.9 


+25.9 
-28.2 


0.98 


±0.25 - 


- ±0 


50 


±0.38 


685.8 


±2.9 


+8.2 
-8.0 


+30.3 
-32.5 


0.98 


±0.50- 


-±o 


75 


±0.63 


779.6 


±3.1 


+12.3 
-12.2 


+32.5 
-34.4 


0.98 


±0.75 - 


-±1 


00 


±0.87 


803.0 


±3.2 


+19.3 
-19.3 


+31.4 
-32.8 


0.99 


±1.00- 


-±1 


25 


±1.12 


784.9 


±3.3 


+19.6 
-19.6 


+29.1 
-30.7 


0.99 


±1.25- 


-±1 


50 


±1.38 


694.9 


±2.9 


+24.0 
-24.0 


+25.1 
-27.0 


1.00 


±1.50- 


-±1 


75 


±1.63 


654.8 


±2.8 


+23.8 
-23.8 


+24.6 
-25.3 


1.00 


±1.75- 


-±2 


00 


±1.88 


592.1 


±2.6 


+16.4 
-16.3 


+22.5 
-23.2 


1.00 


±2.00- 


- ±2 


25 


±2.13 


547.3 


±2.6 


+13.4 
-13.4 


+20.4 
-21.1 


1.00 


±2.25 - 


-±2 


50 


±2.38 


507.2 


±2.4 


+20.2 
-20.2 


+ 18.7 
-19.9 


1.00 



Table 2: The measured differential cross-sections da/dE}f^ and da/drf^^ based 
on the jet algorithm for inclusive-jet photoproduction with —1 < rj^'^^ < 2.5 and 
E^f^ > 17 GeV in the kinematic region given by < 1 GeV^ and 142 < W^p < 
293 GeV. The statistical (5sta.t), uncorrected systematic f^systy' and jet-energy scale 
(Ses) uncertainties are shown separately. The corrections for hadronisation effects 
to be applied to the parton-level NLO QGD calculations (Cy^ad) are shown in the 
last column. 
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„jct 
7 


bill 




\'/ / 


(pb) 


"Stat 


"syst 




"-'had 


-0.75- 


- -0 


50 


-0.56 


4.12 


±0.13 


+0.55 
-0.49 


+0.56 
-0.62 


0.66 


-0.50- 


- -0 


25 


-0.34 


38.12 


±0.60 


+2.51 
-2.45 


+3.47 
-3.88 


0.87 


-0.25 - 


-+0 


00 


-0.11 


108.8 


±1.1 


+3.8 
-3.8 


+7.3 
-7.9 


0.94 


+0.00 - 


-+0 


25 


+0.13 


176.2 


±1.4 


+3.8 

-3.6 


+9.7 
-10.6 


0.97 


+0.25 - 


-+0 


50 


+0.38 


238.3 


±1.7 


+3.3 
-3.2 


+11.8 
-12.7 


0.97 


+0.50 - 


-+0 


75 


+0.63 


282.0 


±1.8 


+4.6 
-4.4 


+13.0 
-14.1 


0.97 


+0.75 - 


-+1 


00 


+0.87 


309.5 


±2.0 


+5.9 
-5.9 


+14.0 
-14.5 


0.98 


+1.00- 


- +1 


25 


+ 1.12 


314.5 


±2.0 


+5.1 
-5.1 


+13.3 
-14.4 


0.99 


+ 1.25- 


- +1 


50 


+ 1.38 


281.6 


±1.8 


+5.9 
-5.9 


+11.5 
-12.2 


1.00 


+1.50- 


- +1 


75 


+ 1.63 


262.9 


±1.7 


+6.5 
-6.5 


+11.2 
-11.4 


1.00 


+ 1.75- 


- +2 


00 


+ 1.88 


225.3 


±1.6 


+4.3 
-4.3 


+9.3 
-9.9 


1.00 


+2.00- 


- +2 


25 


+2.13 


200.9 


±1.5 


+2.6 
-2.6 


+8.4 
-9.1 


1.01 


+2.25 - 


-+2 


50 


+2.38 


173.2 


±1.3 


+5.3 

-5.3 


+7.7 

-8.3 


1.01 



Table 3: The measured differential cross-sections da/drj^^^ based on the kj- jet 
algorithm for inclusive-jet photoproduction with E^rf^ > 21 GeV in the kinematic 
region given by Q'^ < 1 GeV"^ and 142 < W.yp < 293 GeV. Other details as in the 
caption to Table\E 
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E}^' bin 






^stat 


'^syst 




Chad 


(GeV) 


(GeV) 


(pb/GeV) 










-1 < < 


17- 


21 


18.6 


35.99 


±0.17 


-Ul 1 ^ 

-1.14 


-L9 Ozl 

-2.21 


0.93 


21 - 


25 


22.6 


7.522 


±0.071 


+0.337 
-0.331 


+0.543 
-0.592 


0.91 


25 - 


29 


26.5 


1.695 


±0.032 


-0.109 


-t-O 1 ^7 

-0.167 


0.89 


29 - 


35 


30.9 


0.268 


±0.010 


-0.024 


_i_n nQO 

-0.039 


0.85 


35 - 


41 


37.0 


0.0138 


±0.0018 


-f-U.UUZD 

-0.0021 


in 009*^ 

-0.0031 


0.80 








< 










17- 


21 


18.7 


113.61 


±0.31 


1 1 CO 

-hi.uo 

-1.67 


1 A /I c: 
-|-4.40 

-4.72 


0.99 


21 - 


25 


22.7 


37.73 


±0.17 


-0.53 


-1-1 70 

-1.83 


0.98 


25 - 


29 


26.7 


14.27 


±0.11 


-0.24 


-UO 79 

-0.74 


0.98 


29 - 


35 


31.4 


5.034 


±0.052 


-f-u.uyo 
-0.087 


_i_n 9fii 
-0.300 


0.97 


35 - 


41 


37.5 


1.490 


±0.027 


-0.039 


_i_n OQO 
-0.098 


0.96 


41 - 


47 


43.4 


0.485 


±0.015 


+0.016 
-0.015 


+0.031 
-0.032 


0.96 


47 - 


55 


50.1 


0.1356 


±0.0068 


+0.0059 
-0.0056 


+0.0099 
-0.0110 


0.95 


55 - 


71 


59.8 


0.0220 


±0.0017 


in nm A 
-0.0013 


in nm Q 
-0.0019 


0.93 


71 - 


95 


76.5 


0.00075 


±0.00022 


-0.00008 


4-0 0001 

-0.00015 


0.93 








1 < 


r/j^* < 1.5 








17- 


21 


18.7 


55.26 


±0.21 


-2.07 


4-1 R1 

-"1- 1 .Ol 

-1.92 


1.00 


21 - 


25 


22.7 


20.17 


±0.12 


+0.43 
-0.43 


+0.77 
-0.83 


1.00 


25- 


29 


26.7 


8.723 


±0.082 


+0.133 
-0.133 


+0.355 
-0.387 


0.99 


29- 


35 


31.5 


3.461 


±0.044 


+0.062 
-0.060 


+0.173 
-0.184 


0.99 


35- 


41 


37.6 


1.209 


±0.026 


+0.031 
-0.031 


+0.060 
-0.058 


0.97 


41 - 


47 


43.6 


0.486 


±0.016 


+0.013 
-0.013 


+0.026 
-0.027 


0.99 


47- 


55 


50.3 


0.1911 


±0.0089 


+0.0074 
-0.0073 


+0.0117 
-0.0130 


0.97 


55- 


71 


60.5 


0.0454 


±0.0031 


+0.0023 
-0.0023 


+0.0028 
-0.0031 


0.96 


71 - 


95 


77.5 


0.00508 


±0.00083 


+0.00083 
-0.00083 


+0.00041 
-0.00038 


0.99 



Table 4: The measured differential cross-sections da/dE^r^ based on the kr jet 
algorithm for inclusive-jet photoproduction with Ei^^ > 17 GeV in different regions 
j^ct kinematic region given by Q"^ < 1 GeV"^ and 142 < W^p < 293 GeV. 

Other details as in the caption to Table\^ 
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4'* 


bin 






<^stat 


'^syst 


5es 


Chad 


(GeV) 


(GeV) 


(pb/GeV) 










1.5 < r/j^* < 2 


17- 


21 


18.7 


47.51 


±0.19 


+ 1.86 
-1.87 


+ 1.64 
-1.68 


0.99 


21 - 


25 


22.7 


16.70 


±0.11 


+0.48 
-0.48 


+0.66 
-0.68 


1.00 


25 - 


29 


26.7 


6.844 


±0.069 


+0.134 
-0.134 


+0.285 
-0.294 


0.99 


29 - 


35 


31.5 


2.818 


±0.037 


+0.035 
-0.035 


+0.129 
-0.138 


0.99 


35 - 


41 


37.5 


1.055 


±0.023 


+0.020 
-0.020 


+0.056 
-0.052 


0.99 


41 - 


47 


43.7 


0.432 


±0.015 


+0.006 
-0.006 


+0.023 
-0.026 


0.99 


47- 


55 


50.2 


0.1703 


±0.0086 


+0.0020 
-0.0022 


+0.0099 
-0.0119 


0.99 


55 - 


71 


60.6 


0.0334 


±0.0027 


+0.0012 
-0.0011 


+0.0021 
-0.0025 


1.00 


71 - 


95 


78.1 


0.00393 


±0.00081 


+0.00027 
-0.00028 


+0.00041 
-0.00029 


0.98 








2 < 


r/j^* < 2.5 








17- 


21 


18.6 


42.88 


±0.18 


+1.63 

-1.63 


+1.40 
-1.43 


1.00 


21 - 


25 


22.7 


13.479 


±0.097 


+0.325 
-0.325 


+0.523 
-0.591 


1.01 


25 - 


29 


26.7 


5.223 


±0.058 


+0.134 
-0.134 


+0.240 
-0.222 


1.01 


29 - 


35 


31.5 


1.977 


±0.029 


+0.027 
-0.027 


+0.091 
-0.114 


1.00 


35 - 


41 


37.6 


0.708 


±0.018 


+0.016 
-0.015 


+0.039 
-0.038 


1.00 


41 - 


47 


43.6 


0.268 


±0.011 


+0.003 
-0.004 


+0.017 
-0.017 


1.00 


47- 


55 


50.2 


0.0928 


±0.0059 


+0.0065 
-0.0065 


+0.0061 
-0.0082 


1.00 


55 - 


71 


60.3 


0.0192 


±0.0019 


+0.0007 
-0.0007 


+0.0016 
-0.0013 


0.97 


71 - 


95 


74.9 


0.00238 


±0.00065 


+0.00023 
-0.00023 


+0.00024 
-0.00021 


1.11 



Table 5: Continuation of Table\^ 
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E^r^^ bin 






<^stat 


^syst 


^ES 


Chad 


(GeV) 


(GeV) 


(pb/GeV) 










anti-Zcy 


17- 


21 


18.6 


276.84 


±0.50 


+5.55 
-5.55 


+ 11.09 
-11.75 


0.93 


21 - 


25 


22.7 


89.51 


±0.28 


+1.36 
-1.35 


+4.03 
-4.30 


0.93 


25- 


29 


26.7 


34.54 


±0.18 


+0.49 
-0.49 


+1.68 
-1.76 


0.94 


29- 


35 


31.4 


12.771 


±0.091 


+0.180 
-0.170 


+0.640 
-0.714 


0.94 


35- 


41 


37.5 


4.218 


±0.052 


+0.090 
-0.087 


+0.235 
-0.243 


0.94 


41 - 


47 


43.5 


1.567 


±0.032 


+0.030 
-0.028 


+0.093 
-0.101 


0.95 


47- 


55 


50.2 


0.550 


±0.017 


+0.018 
-0.018 


+0.037 
-0.039 


0.95 


55- 


71 


60.3 


0.1139 


±0.0058 


+0.0039 
-0.0034 


+0.0074 
-0.0081 


0.93 


71 - 


95 


77.5 


0.0105 


±0.0014 


+0.0006 
-0.0007 


+0.0009 
-0.0010 


0.96 








SIS 


cone 








17- 


21 


18.7 


278.01 


±0.51 


+4.31 
-4.30 


+10.89 
-11.52 


0.80 


21 - 


25 


22.7 


90.82 


±0.28 


+1.22 
-1.21 


+3.96 
-4.27 


0.81 


25- 


29 


26.7 


35.27 


±0.18 


+0.49 
-0.48 


+1.65 
-1.75 


0.82 


29- 


35 


31.4 


13.059 


±0.090 


+0.185 
-0.179 


+0.651 
-0.695 


0.84 


35- 


41 


37.5 


4.330 


±0.051 


+0.108 
-0.107 


+0.233 
-0.252 


0.85 


41 - 


47 


43.6 


1.639 


±0.031 


+0.031 
-0.029 


+0.090 
-0.100 


0.86 


47- 


55 


50.3 


0.565 


±0.016 


+0.020 
-0.020 


+0.034 
-0.039 


0.87 


55- 


71 


60.5 


0.1199 


±0.0055 


+0.0045 
-0.0042 


+0.0085 
-0.0085 


0.86 


71 - 


95 


77.7 


0.0108 


±0.0013 


+0.0008 
-0.0008 


+0.0012 
-0.0013 


0.90 



Table 6: The measured differential cross-sections da/dE^r^ based on different jet 
algorithms for inclusive-jet photoproduction with —1 < rf^^ < 2.5 in the kinematic 
region given by Q"^ < 1 GeV"^ and 142 < W^p < 293 GeV. Other details as in the 
caption to Table\^ 
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<^stat 




^ES 


Chad 








(pb) 












-0.75 - 


-0.50 


-0.58 


65.63 


±0.84 


+4.16 

-3.88 


+5.83 

-6.51 


0.74 


-0.50- 


-0.25 


-0.35 


200.3 


±1.6 


+8.1 
-7.9 


+13.3 
-14.1 


0.87 


-0.25 - 


+0.00 


-0.12 


369.5 


±2.1 


+10.1 
-9.9 


+20.3 
-21.4 


0.92 


+0.00 - 


+0.25 


+0.13 


518.6 


±2.5 


+7.5 
-7.3 


+25.5 
-26.3 


0.93 


+0.25- 


+0.50 


+0.38 


634.6 


±2.8 


+6.9 
-6.6 


+28.3 
-30.3 


0.93 


+0.50 - 


+0.75 


+0.63 


722.2 


±3.0 


+9.6 
-9.5 


+30.8 
-32.7 


0.93 


+0.75 - 


+ 1.00 


+0.87 


748.2 


±3.1 


+15.5 
-15.5 


+30.0 
-32.1 


0.94 


+1.00- 


+ 1.25 


+ 1.12 


732.3 


±3.2 


+18.3 
-18.3 


+28.4 
-29.8 


0.95 


+1.25- 


+ 1.50 


+ 1.38 


649.7 


±2.8 


+19.1 
-19.1 


+24.0 
-25.8 


0.94 


+1.50- 


+ 1.75 


+1.63 


610.9 


±2.7 


+20.1 
-20.2 


+23.4 
-25.0 


0.94 


+1.75- 


+2.00 


+ 1.87 


569.5 


±2.6 


+16.2 
-16.2 


+22.4 
-24.1 


0.94 


+2.00 - 


+2.25 


+2.13 


536.2 


±2.6 


+16.1 
-16.1 


+21.3 
-22.8 


0.94 


+2.25 - 


+2.50 


+2.38 


488.6 


±2.5 


+18.5 
-18.5 


+18.7 
-21.2 


0.94 


SIScone 


-0.75- 


-0.50 


-0.58 


69.65 


±0.87 


+4.67 
-4.36 


+6.21 
-7.06 


0.66 


-0.50- 


-0.25 


-0.35 


205.3 


±1.6 


+8.7 
-8.7 


+13.5 
-14.7 


0.77 


-0.25 - 


+0.00 


-0.12 


379.3 


±2.2 


+11.0 
-10.9 


+20.4 
-22.1 


0.81 


+0.00 - 


+0.25 


+0.13 


527.7 


±2.6 


+8.3 
-8.0 


+25.6 
-27.1 


0.82 


+0.25- 


+0.50 


+0.38 


645.1 


±2.9 


+6.8 
-6.7 


+28.1 
-30.6 


0.82 


+0.50 - 


+0.75 


+0.63 


731.6 


±3.1 


+9.1 
-9.0 


+30.1 
-32.5 


0.82 


+0.75 - 


+1.00 


+0.87 


756.8 


±3.2 


+13.8 
-13.8 


+29.9 
-31.8 


0.83 


+1.00- 


+ 1.25 


+ 1.12 


735.0 


±3.3 


+14.7 
-14.6 


+27.7 
-29.3 


0.83 


+ 1.25- 


+ 1.50 


+ 1.38 


650.0 


±2.9 


+15.2 
-15.3 


+23.3 
-24.3 


0.82 


+ 1.50- 


+ 1.75 


+ 1.63 


598.7 


±2.7 


+14.3 
-14.3 


+22.3 
-23.1 


0.81 


+ 1.75- 


+2.00 


+ 1.87 


549.0 


±2.6 


+10.8 
-10.9 


+20.9 
-22.2 


0.81 


+2.00 - 


+2.25 


+2.13 


529.8 


±2.5 


+14.1 
-14.1 


+20.5 
-22.0 


0.80 


+2.25 - 


+2.50 


+2.38 


518.5 


±2.5 


+24.3 

-24.4 


+20.1 
-21.6 


0.79 



Table 7: The measured differential cross-sections da/drj^'^^ based on different jet 
algorithms for inclusive-jet photoproduction with E^f^ > 17 GeV in the kinematic 
region given by Q"^ < 1 GeV"^ and 142 < Wyp < 293 GeV. Other details as in the 
caption to Table\E 
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anti-/cr 


SIScone 


experimental uncertainties 


jet energy scale 


+1.8cr/ 
-1.7/0 


+1.80/ 
-1.8/0 


+1.70/ 
-1.6 /o 


luminosity 


±0.6% 


±0.6% 


±0.6% 


uncorrelated 


+0.3C/ 
-0.4 /O 


+0.30/ 
-0.4/0 


+0.30/ 
-0.4 /O 


statistical 


±0.2% 


±0.2% 


±0.2% 


theoretical uncertainties 


terms beyond NLO 


+2.40/ 


+2.30/ 
-2.4/0 




photon PDFs 


+2.30/ 
-0.9/0 


+2.20/ 
-0.9/0 




proton PDFs 


±1.0% 


±1.0% 


±1.0% 


hadronisation 


±0.4% 


±0.4% 


±0.2% 


total uncertainty 




+4.00/ 
-3.4/0 


+3.90/ 
-3.4/0 


+4.20/ 
-4.0/0 



Table 8: Experimental, theoretical and total uncertainties in the determination 
of as{Mz) from the kx, anti-kx and SIScone analyses. 







'^uncorr 


'^corr 




(GeV) 










22.7 


0.1561 


±0.0011 


+0.0035 
-0.0048 


+0.0106 
-0.0089 


26.7 


0.1493 


±0.0007 


+0.0033 
-0.0034 


+0.0083 
-0.0070 


31.4 


0.1443 


±0.0005 


+0.0035 
-0.0030 


+0.0069 
-0.0059 


37.5 


0.1396 


±0.0007 


+0.0032 
-0.0031 


+0.0057 
-0.0051 


43.6 


0.1359 


±0.0011 


+0.0032 
-0.0030 


+0.0051 
-0.0047 


50.2 


0.1328 


±0.0014 


+0.0037 
-0.0034 


+0.0047 
-0.0045 


60.3 


0.1283 


±0.0024 


+0.0040 
-0.0036 


+0.0041 
-0.0041 



Table 9: The ag values determined in each (-E^*) value from the analysis of the 
measured da/dE^^^ cross section based on the kx jet algorithm. The uncorrelated 
r^uncorry' O'T-C? Correlated f^corrj experimental and theoretical (6th) uncertainties are 
listed separately. 
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Figure 1: Overview of the relative theoretical uncertainties for inclusive-jet cross 
sections in photoproduction in the kinematic region of the measurements as func- 
tions of (a,h) E}^^ and (c,d) rf^^ for the kx, anti-k^ and SIScone jet algorithms. 
Shown are the relative uncertainties induced by the terms beyond NLO, the proton 
PDFs, the value of as{Mz), the modelling of the QCD cascade and hadronisation 
and the photon PDFs. 
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Figure 2: (a) The measured differential cross-section da/dE^f^ based on the 
jet algorithm for inclusive-jet photoproduction with —1 < rj^^^ < 2.5 (dots) in the 
kinematic region given by Q"^ < 1 GeV"^ and 142 < W^p < 293 GeV. The NLO 
QCD calculation (solid line), corrected to include hadronisation effects and using 
the ZEUS-S (GRV-HO) parameterisations of the proton (photon) PDFs, is also 
shown, (b) The relative difference between the measured da/dE}^^ and the NLO 
QCD calculation (dots). In both figures, the inner error bars represent the statistical 
uncertainties; the outer error bars show the statistical and systematic uncertainties 
not associated with the uncertainty in the absolute energy scale of the jets, added 
in quadrature; the shaded band displays the uncertainty due to the absolute energy 
scale of the jets and the hatched band displays the total theoretical uncertainty. In 
some bins, the error bars on the data points are smaller than the marker size and 
are therefore not visible. 
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Figure 3: (a) The measured differential cross-section da/drf'^^ based on the kx 
jet algorithm for inclusive-jet photoproduction with E^rp" > 17 GeV (dots) in the 
kinematic region given by Q"^ < 1 GeV"^ and 142 < W^p < 293 GeV. (b) The 
relative difference between the measured da/drj^'^^ and the NLO QGD calculation 
(dots). Other details as in the caption to Fig.\^ 
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Figure 4: The measured differential cross-sections (a) da/dE^^^ and (h) da/drj^^^ 
based on the jet algorithm for inclusive-jet photoproduction with E^r^^ > 17 GeV 
and —1 < r/j'^* < 2.5 (dots) in the kinematic region given by < 1 GeV"^ and 
142 < Wryp < 293 GeV. For comparison, the NLO QGD calculations including an 
estimation of non-perturbative effects (see text) are also shown. Other details as in 
the caption to Fig. [H 
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Figure 5: The measured differential cross-section da/drj^'^^ based on the jet 
algorithm for inclusive-jet photoproduction with E-!f'^ > 21 GeV (dots) in the kine- 
matic region given by Q'^ < I GeV'^ and 142 < W^p < 293 GeV. In (h), the NLO 
QCD calculations including an estimation of non-perturbative effects (see text) are 
also shown. Other details as in the caption to Fig. [3 
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Figure 6: The measured differential cross-sections (a) da/dE^^^ and (h) da/drj^^^ 
based on the kx jet algorithm for inclusive-jet photoproduction with E^r^^ > 17 GeV 
and —1 < r/j'^* < 2.5 (dots) in the kinematic region given by < 1 GeV^ and 
142 < Wjp < 293 GeV. For comparison, the NLO QGD calculations using different 
parameterisations of the photon PDFs are also shown. Other details as in the 
caption to Fig. [3 
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Figure 7: The measured differential cross-sections (a) da/dE^f^ and (h) da/dr]^'^^ 
based on the kr jet algorithm for inclusive-jet photoproduction with E^rf^ > 17 GeV 
and —1 < rf^^ < 2.5 (dots) in the kinematic region given by Q"^ < 1 GeV"^ and 
142 < W^p < 293 GeV. For comparison, the NLO QGD calculations using different 
parameterisations of the proton PDFs are also shown. Other details as in the 
caption to Fig. 
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Figure 8: The measured differential cross-sections da/dE^f^ based on the kr jet 
algorithm for inclusive-jet photoproduction with E^rf^ > 17 GeV in different regions 
ofrf'^^ (dots) in the kinematic region given by Q"^ < 1 GeV^ and 142 < W^p < 293 
GeV. Each cross section has been multiplied by the scale factor indicated in brackets 
to aid visibility. Other details as in the caption to Fig. O 
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Figure 9: The relative differences between the measured differential cross-sections 
da/dE^r^^ presented in Fig. and the NLO QCD calculations (dots). Other details 
as in the caption to Fig. 
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Figure 10: The relative differences between the measured differential cross-sections 
da/dE}^^ presented in Fig. and the NLO QCD calculations (dots). The relative 
differences between the predictions based on the calculations including an estimation 
of non-perturbative effects (see text) and the NLO QCD calculation are also shown. 
Other details as in the caption to Fig. [B 
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Figure 11: The relative differences between the measured differential cross-sections 
da/dE^f^ presented in Fig. and the NLO QCD calculations (dots). The relative 
differences between the predictions based on different photon PDFs and that based 
on the ZEUS-S/GRV-HO sets are also shown. Other details as in the caption to 
Fig.lE 
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Figure 12: The relative differences between the measured differential cross-sections 
da/dE^f'^ presented in Fig. and the NLO QCD calculations (dots). The relative 
differences between the predictions based on different proton PDFs and that based 
on the ZEUS-S/GRV-HO sets are also shown. Other details as in the caption to 
Fig.\B 
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Figure 13: The measured differential cross-sections (a) da/dE^r^ and (h) 
da/drj^'^^ based on different jet algorithms for inclusive-jet photoproduction with 
E^^ > 17 GeV and -1 < 77-''^* < 2.5 (dots) in the kinematic region given by 
< 1 GeV"^ and 142 < W^p < 293 GeV. The anti-kx and kx cross sections 
were multiplied by the scale factors indicated in brackets to aid visibility. The lower 
part of the figure shows the hadronisation correction factors applied to the NLO 
calculations together with their uncertainty (hatched bands) for each jet algorithm; 
the hadronisation correction factor for the kx algorithm was shifted by the value 
indicated in brackets to aid visibility. Other details as in the caption to Fig. 
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Figure 14: T/ie ratios between the measured cross-sections (a) da/dE^^^ and (h) 
dojdri^^^ and the NLO QCD calculations (dots) from Fig. [73 Other details as in 
the caption to Fig. [H 
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Figure 15: The ratios of the measured cross-sections anti-kr/kT, SIScone/kT 
and anti-k-T /SIScone (dots) as functions of (a) E}^^ and (b) rj^^^. In these plots, 
the outer error bars also include the uncertainty on the absolute energy scale of the 
jets. The predicted ratios based on calculations which include up to O(a^) terms 
are also shown (solid lines). The hatched bands display the theoretical uncertainty 
on the ratio. The dashed lines indicate the ratios of the hadronisation correction 
factors and the dash-dotted lines represent the ratios of the NLO QCD calculations 
including an estimation of non-perturbative effects (see text). Other details as in 
the caption to Fig. 
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Jet production at high Q 
HI (Eur Phys J C 65 (2010) 363) 
Jet production at low Q 
HI (Eur Phys J C 67 (2010) 1) 

Inclusive-jet cross sections in pp 
CDF (Phys Rev Lett 88 (2002) 042001) 

Inclusive-jet cross sections in pp 
DO (Phys Rev D 80 (2009) 111107) 

Event shapes in e^e" at NNLO+NLLA 
OPAL (Eur Phys J C 71 (2011) 1733) 

Event shapes in e^e" 
L3 (Phys Lett B 536 (2002) 217) 

Jet rates in e^e" 

DELPHI (Eur Phys J C 38 (2005) 413) 

HERA average 2004 
(C. Glasman, hep-ex/0506035) 

HERA combination 2007 
(Hlprelim-07-132/ZEUS-prel-07-025) 

World average 2009 
(S. Bethke, Eur Phys J C 64 (2009) 689) 
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Figure 16: Extracted as{Mz) value from this analysis (upper dot). For com- 
parison, determinations from other experiments and reactions, the HERA average 
2004, the HERA combination 2007 and the world average 2009 are also shown. The 
horizontal error bars represent the experimental and theoretical uncertainties added 
in quadrature. The shaded hand represents the uncertainty of the world average. 
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Figure 17: The as values determined in each (i?^*) value from the analysis of the 
measured da/dE^r^^ cross section based on the kr jet algorithm (open circles). The 
error bars represent the uncorrelated experimental uncertainties; the shaded area 
represents the correlated experimental uncertainties and the hatched area represents 
the correlated experimental and theoretical uncertainties added in quadrature. The 
solid line indicates the renormalisation- group prediction at two loops obtained from 
the corresponding as{Mz) value determined in this analysis. 
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